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Abstract. We review theoretical models of Population III. 1 star formation, focusing on the proto- 
stellar feedback processes that are expected to terminate accretion and thus set the mass of these 
stars. We discuss how dark matter annihilation may modify this standard feedback scenario. Then, 
under the assumption that dark matter annihilation is unimportant, we predict the mass of stars form- 
ing in 12 cosmological minihalos produced in independent numerical simulations. This allows us to 
make a simple estimate of the Pop III. 1 initial mass function and how it may evolve with redshift. 
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INTRODUCTION 

The first stars laid the foundations for galaxy formation. As described below, we expect 
that these Population III. 1 stars, i.e. having negligible metallicity and experiencing 
negligible influence from other astrophysical sources external to their minihalos, were 
massive and thus injected significant radiative, mechanical and chemical feedback into 
their surroundings. We need to understand the processes that set the mass of these stars 
in order to have a complete theory of galaxy formation and evolution. Supermassive 
black holes reside in the centers of most large galaxies and some massive examples have 
been seen at very early times [1, 2]. The remnants of the first stars may have been the 
seeds of these black holes. 

We first review the accretion and feedback processes that are expected to occur during 
the birth of Pop III. 1 stars. We emphasize some of the uncertainties in existing models. 
We then describe preliminary work applying these models to a population of Pop III. 1 
minihalos in order to predict the initial mass function (IMF) of the stars they produce. 



ACCRETION AND FEEDBACK PROCESSES 

Here we give a brief outline of these processes, which have been reviewed in more detail 
by Tan & McKee[3]. Small dark matter halos grow and merge together to eventually 
form ~ 1O 6 M ' 'mini-halos" by redshifts z ~ 30 — 20. Baryons collect in these halos 
and cool down to ~ 200 K due to the presence of trace amounts of H2, the formation 
of which is catalyzed by small concentrations of free electrons via the H~ channel. 
The baryons start to dominate the mass density in the innermost ~ 1 parsec. Above 



the critical density of «h ~ 10 4 cm -3 , cooling is inefficient, contraction is slow and 
subsonic, and the temperature gradually rises towards the center. The gas develops a self- 
similar, power law density structure with p oc r~ 2 - 2 , i.e. slightly steeper than a singular, 
isothermal sphere. When the central density exceeds ~ 10 10 cm~ 3 , the gas becomes 
fully molecular via the 3-body formation channel, cooling is much more efficient and 
dynamical collapse is initiated from the inside-out. 

An analytic estimate of the accretion rate to the central star+disk system is [4]: 

m %d = O.O26e^ ,15 / 7 (M/M )- 3 / 7 M Q yr~\ (1) 

where M is the mass originally enclosed by the currently accreting mass shell, is the 
instantaneous star formation efficiency, and K' is the "entropy parameter" normalizing 
the polytropic description of the initial core, P = Kp Y with y ~ 1.1, and is defined as 

P/P J _ T _k f l0 4 cm- 3 
1.88xl0 12 cgs 300 KV n H 

where T' s is an effective temperature that includes the modest effect of subsonic turbu- 
lent motions that are seen in the numerical simulations of Abel, Bryan & Norman[5]. 

The gas core has significant rotational motions [5] and so a disk should form inside the 
sonic point of the accretion flow. We expect the disk mass to build up until gravitational 
instabilities produce spiral density waves and other structures that allow transfer of mass 
inwards and angular momentum outwards through the disk. For fiducial accretion rates 
and disk sizes, we expect the optically thick region of the disk to be stable to gravitational 
fragmentation^]. Most simulations of Pop III. 1 star formation appear to produce single 
stars [7], although these have followed only the earliest stages of accretion when the 
central protostar has yet to accrete most of its expected final mass. 

The accretion disk is the first place where turbulent eddies have time to amplify 
significantly seed magnetic field produced by the Biermann battery mechanism during 
the formation of the halo[6]. If such an accretion disk dynamo produces large-scale, 
dynamically- strong B-fields, then a bipolar magneto-centifugally-driven outflow will be 
launched, which can reduce e^. Tan & Blackman[6] estimated ~ 0.5 by the time 
m* ~ 200M Q , which is not as important as radiative feedback in the fiducial case with 
K' = 1 (see below). 

Matter joins the protostar via an accretion boundary layer shock and the internal 
energy of the gas that is advected into the protostar is set by the properties of this shock, 
the accretion rate and the amount of energy that is radiated. The internal energy of the 
accreted gas in turn helps set the size of the protostar[8, 4]. 

Weakly interacting massive particle (WIMP) dark matter may self-annihilate and 
provide significant energy input into the protostar[9, 10], stabilizing its initial radius 
at ~ 10 AU, rather than ~ 0.1 AU in the fiducial case. The subsequent evolution is 
uncertain[ll], but contraction to the main sequence may be delayed with respect to 
the fiducial case without such support (see below), affecting the mass scale at which 
radiative feedback truncates accretion. Uncertainties include the assumed central density 
profile of dark matter and how this may be affected by supersonic baryonic infall and 
disk accretion, including heating by spiral density waves. For the remainder of this 
article, we consider the case of negligible WIMP annihilation heating. 




The protostellar evolution can be followed given knowledge of the accretion rate[4, 
12]. Once the protostar is older than its local Kelvin-Helmholz time, it starts contracting 
towards the main sequence, typically settling there by ~ 1OOM . The star could continue 
to accrete in this state if feedback is not strong enough to reverse accretion. Given the 
accretion rates of equation (1) and stellar evolution timescales of Schaerer[13], the max- 
imum accreted stellar mass before supernova explosion is m* 5max — 19OOi^ /L28 M [14]. 

McKee & Tan[14] considered a number of feedback processes that will intervene 
before a supernova occurs. Destruction of H2 coolant in the minihalo by FUV photons 
was not found to reduce the accretion rate to the star significantly, since by the time 
this occurs the protostar has already developed a potential well deep enough for atomic 
cooling to allow accretion. Lyman- a radiation pressure was also found to be ineffective 
since it requires multiple scatterings to build up dynamically important pressures, but 
this leads to diffusion and escape of the photons along bipolar cavities. 

Direct ionization of the accretion flow to create a ~ 25, 000 K HII region whose high 
thermal pressure can reverse accretion is the first important feedback process, occurring 
by ~50 Mq in the fiducial case. The HII region quickly breaks out in all directions 
above and below the disk, reducing by factors of several. By the time m* = 1OOM 
the accretion rate from disk-shielded directions is about a factor of 10 smaller than in 
the no-feedback case. These estimates use a ID model for the disk vertical thickness and 
assume gas from disk-shielded directions accretes as it would in the absence of feedback. 

The protostar ionizes the accretion disk atmosphere, which redirects some ionizing 
flux to outer parts of the disk that are neutral. An ionized, "photoevaporative" flow 
develops from these outer regions with mass-loss rate[15] 

m evap ~ 4.1 x \Q- 5 S% 2 T$ A m l Jl 2 M Q yr" 1 , (3) 

where S49 is the H-ionizing photon luminosity in units of 10 49 photons s _1 , T4 is the 
ionized gas temperature in units of 10 4 K, and the star+disk mass has been normalized 
to 1OOM . This analytic estimate has not yet been confirmed by numerical simulation. 
The final accreted mass, m*y, is estimated by the condition m evap = m* d , yielding[14]: 

m* f = 145^ 6 °/ 47 (2.5/r 4 ) o - 24 (/ sh /O.2) 28 / 47 (e w /O.25) 12 / 47 M , (4) 

where / s h is the shadowing factor, i.e. fraction of the sky seen from the protostar 
blocked by the disk, and is the mass-averaged value of over the formation of 
the star. The fiducial normalizations of these parameters were derived from calculations 
of the protostellar and accretion disk evolution coupled with the effects of ionizing 
feedback. The predicted mass of Pop III. 1 stars in the fiducial case is ~ 15OM . The 
strongest parameter dependence is with K', i.e. via the accretion rate. High accretion 
rate protostars reach higher masses because photoevaporative truncation occurs at a 
higher mass, partly because contraction to the main sequence and HII region breakout 
were delayed, thus raising e^, and partly because the star must produce more ionizing 
photons to increase photoevaporative mass-loss to counter the higher accretion rate. 

A number of factors are likely to raise this estimate of the final accreted mass, 
including instabilities during the breakout of the HII region and rotation of the protostar 
that reduces the equatorial flux of ionizing photons and thus the photoevaporative mass- 
loss. 



THE POP IIL1 IMF 



O'Shea & Norman [16] studied the properties of Pop III. 1 pre-stellar cores as a function 
of redshift. They found that cores at higher redshift are hotter in their outer regions, have 
higher free electron fractions and so form larger amounts of H2 (via H ), although these 
are always small fractions of the total mass. As the centers of the cores contract above 
the critical density of 10 4 cm 3 , those with higher H2 fractions are able to cool more 
effectively and thus maintain lower temperatures to the point of protostar formation. The 
protostar thus accretes from lower-temperature gas and the accretion rates, proportional 
to c 3 °c T 3 / 2 , are smaller. Measuring infall rates at the time of protostar formation at 
the scale of M = 1OOM , O'Shea & Norman found accretion rates of ~ 1O~ 4 M yr _1 
at z = 30, rising to ~ 2 x 1O~ 2 M yr _1 at z = 20. They used these accretion rates to 
estimate the mass of the star that forms by finding the condition of when the accretion 
timescale was equal to the Kelvin-Helmholtz timescale, thus predicting an increase in 
the typical Pop III. 1 stellar mass as z decreases over this range. 

Here we present preliminary results of a more detailed estimate of the Pop III. 1 
stellar masses that result from these cores. First we estimate the expected protostellar 
accretion rates in the O'Shea & Norman[16] cores (Fig. 1). In the simulations, the inner 
regions have not yet reached the densities associated with dynamical collapse, so the 
observed mass infall rate underestimates the expected accretion rates. Thus, we assume 
the solution describing the collapse of a singular polytropic sphere (SPS)[4] holds for 
the inner mass shells. We derive this by evaluating the interior mass averaged K' profiles 
of the simulated cores and then using eq. 1. The infall mass flux through spherical shells 
provides another method of estimating the eventual accretion rate to the protostar. For 
Hunter's mildly subsonic solution [17] used in eq. 1 the accretion rate increases from 
0.70c 3 /G at large radii (r ^> c s t, where t = is the moment of protostar formation) to 
2.58c 3 . /G at small radii, an increase of a factor 3.7. Thus we increase the observed infall 
mass flux by this factor, and refer to this as the "enhanced infall rate". Beyond the mass 
shell where this mass accretion rate becomes larger than that predicted by the collapse of 
the SPS, we assume the accretion rate is equal to their geometric mean. This is designed 
to approximate the transition from the SPS solution to the infall solution. In the outer 
parts of the core, the gas has not yet virialized, K' becomes large and the SPS-based 
accretion rate becomes unrealistically large, exceeding the enhanced infall rate. In this 
case, we adopt only the enhanced infall rate as the estimate of the mass accretion rate. 

We then estimate the final protostellar mass, m*j, using the feedback models of 
McKee & Tan[14]. The collapsed mass at HII region breakout can be evaluated from 
their eq. (41) with parameters: 0E dd = 0.3, /i = l/y/l, f d = 1/3, 0s = 1, T4 = 2.5, 

/ Kep = 0.5, mii = 2r g and £, d = 1, then m, ;break = 48AK' 25 / l4 M Q = 6A2ml d /2n 3 M Q . 
Following the accretion history of each halo, we assume that once this condition is met 
HII region breakout is fast and m* d is immediately reduced by the disk shadowing factor 
/ S h = 0.2. Finally, we calculate when disk photoevaporation overwhelms accretion, i.e. 

^evap — > 1-96 x IQr^mJ^ M©yr _1 > m w /[14], at which point we assume the protostar 
has reached its m^f. A more detailed analysis, fully coupling m* d to these protostellar 
evolution and feedback models will be presented in a future paper. 

Figure 2 shows the derived final protostellar masses as a function of the collapse 




FIGURE 1. Estimates of mass accretion rates of the 12 Pop III. 1 cores of O'Shea & Norman[16], each 
indicated by a different line style and color, with fiducial rate of Tan & McKee[4] (eq. 1) shown by the 
black dot-dashed line: (a) Top left: Accretion rate based on SPS collapse with K' profiles (averaged over 
interior mass) (eq. 1); (b) Top right: Accretion rate based on mass flux through spherical shells at the end 
of the simulations; (c) Bottom left: Accretion rate as in (b), now smoothing at 50% weight with the two 
adjacent mass shells and an enhancement by a factor of 3.7 as expected if the Hunter[17] solution holds; 
(d) Bottom right: Final adopted accretion rate, assuming inner SPS solution (a), then geometric mean of 
the (a) and (c) results when the latter starts to dominate, then (c) if (a) starts to exceed (c) again (see text). 



redshift of the minihalo. There is a weak trend for lower-mass stars to be formed at 
higher redshift due to lower accretion rates, as suggested by O'Shea & Norman[16]. 
However, there is significant dispersion about this trend, since although the high z 
protostars do initially have low accretion rates, these can increase due to large mass infall 
rates at the mass shells for M ~ 100 — 1000 M & . Figure 2 also shows the distribution of 
the initial Pop III. 1 stellar masses, i.e. their IMF, based on these 12 minihalos. These 
are the first halos to collapse in their local cosmological volumes, so may be a biased 
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FIGURE 2. Left: Final protostellar mass due to photoevaporative feedback, m*/, versus z for the 12 
halos of O'Shea & Norman[16]. Right: IMF of the stars forming in these Pop III. 1 minihalos. 



sample. The sample is also small, so we note only some basic features: we see very 
few stars below 80 M , a steep rise in the IMF around 100 M , then a general decline to 
higher masses, with the most massive star being ~ 600 M . The mean mass is ~ 25OM , 
the median is ~ 215 M Q and the dispersion is ~ 16OM . 
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